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Abstract 

This paper explores of the role of unitary braiding operators in quantum com- 
puting. We show that a single specific solution R (the Bell basis change 
matrix) of the Yang-Baxter Equation is a universal gate for quantum com- 
puting, in the presence of local unitary transformations. We show that this 
same R generates a new non-trivial invariant of braids, knots, and links. 
Other solutions of the Yang-Baxter Equation are also shown to be universal 
for quantum computation. The paper discusses these results in the context of 
comparing quantum and topological points of view. In particular, we discuss 
quantum computation of link invariants, the relationship between quantum 



entanglement and topological entanglement, and the structure of braiding in 
a topological quantum field theory. 

1 Introduction 

It is a challenge to unravel the relationships among quantum entanglement, 
topological entanglement and quantum computation. In this paper, we show 
some of the pieces in this puzzle and how they fit together. In no way do 
we claim to have assembled the entire puzzle! That is a challenge for futher 
work. In order to introduce our problems, and explain what we have done 
with them, the next few paragraphs will give capsule summaries of each of 
the major points of view taken in this study. We then describe in more 
detail what is contained in each separate section of the paper. The paper 
itself strives to be self-contained, and to describe carefully the issues involved, 
particularly with topological structures that may be unfamiliar to a physics 
audience. 

Quantum computing can be regarded as a study of the structure of the 
preparation, evolution and measurement of quantum systems. In the quan- 
tum computation model, an evolution is a composition of unitary transforma- 
tions (finite dimensional over the complex numbers). The unitary transfor- 
mations are applied to an initial state vector that has been prepared for this 
process. Measurements are projections to elements of an orthonormal basis 
of the space upon which the evolution is apphed. The result of measuring a 
state l'^), written in the given basis, is probabilistic. The probability of ob- 
taining a given basis element from the measurement is equal to the absolute 
square of the coefficient of that basis element in the state being measured. 

It is remarkable that the above lines constitute an essential summary of 
quantum theory. All applications of quantum theory involve filling in details 
of unitary evolutions and specifics of preparations and measurements. 

One hopes to build powerful quantum computers. Such hopes would be 
realized if there were reliable ways to implement predetermined patterns of 
unitary evolution and measurement. In the course of trying to understand the 
potential for quantum computing, it became apparent that arbitrary finite 
dimensional unitary transformations can be built from a relatively small set 



2 



of primitives. A standard set of primitives consists in all two-dimensional 
unitary transformations, together with a choice of one sufficiently robust 
four-dimensional transformation such as the CNOT gate discussed in the 
first section of this paper. One says that CNOT, together with single qubit 
gates (two dimensional unitary transformations) is universal for quantum 
computation. 

Probability in quantum mechanics acts quite differently than classical 
probability. Entangled quantum states embody this difference. An example 
of an entangled state is the two-qubit state \'4>) = (|00) + \ll))/\/2. This 
state is not decomposable as a tensor product of single-qubit states, and 
a measurement in one of its tensor factors will determine the outcome in 
the other factor. Implicit in entanglement is the phenomenon of quantum 
non-locality: physical access to the measurement of one tensor factor or 
the other may be separated by an abitrary spatial interval. The result of a 
measurement can have the appearance of instantaneous determination across 
an arbitrary distance. 

Entanglement and quantum computing are related in a myriad of ways, 
not the least of which is the fact that one can replace the CNOT gate 
by another gate R and maintain universality (as described above) just so 
long as R can entangle quantum states. That is, R can be applied to some 
unentangled state to produce an entangled state. It is of interest to examine 
other sets of universal primitives that are obtained by replacing CNOT by 
such an R. 

Contemplating the inherent non-locality of entangled states, it is natural 
to ask whether there are relationships between topological entanglement and 
quantum entanglement. Topology studies global relationships in spaces, and 
how one space can be placed within another, such as knotting and linking 
of curves in three-dimensional space. One way to study topological entan- 
glement and quantum entanglement is to try making direct correspondences 
between patterns of topological linking and entangled quantum states. One 
approach of this kind was initiated by Aravind as we discuss in section 8 of 
this paper and also in [26, 27]. A deeper method (we believe) is to consider 
unitary gates R that arc both universal for quantum computation and are 
also solutions to the condition for topological braiding. Such matrices R are 
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unitary solutions to the Yang-Baxter equation, as explained in section 2. We 
are then in a position to compare the topological and quantum properties of 
these transformations. In this way, we can explore the apparently complex 
relationship among topological entanglement, quantum entanglement, and 
quantum computational universality. It is this exploration that is the theme 
of this paper. 

In this paper, we prove that certain solutions to the Yang-Baxter equation 
together with local unitary two dimensional operators form a universal set of 
quantum gates. In the first version of this result, we generate CNOT using 
a solution to the algebraic Yang-Baxter equation. In the second version, 
we generate CNOT using versions of the braiding Yang-Baxter equation. 
Results of this kind follow from general results of the Bryhnskis [6] about 
universal quantum gates. Here, we give explicit proofs by expressing the 
CNOT gate in terms of solutions to the Yang-Baxter equation (and local 
unitary transformations). 

Section 2 of the paper defines the Yang-Baxter equation, gives unitary 
examples and proves the results about universal gates. We regard these 
results as a significant elementary step in relating quantum topology and 
quantum computing. The results say that quantum computing can be framed 
in the context of quantum topology. They also say that quantum computing 
can be framed in those statistical mechanics contexts where the solutions to 
the Yang-Baxter equation are natural structures. 

Certainly the Yang-Baxter Equation is a natural structure in thinking 

about the topology of braids, knots, and links. In section 3, we formalize an 
extension of the Artin Braid group that can accomodate local operators so 
that this extended braid group can represent any unitary transformantion. 
Section 3 shows how to use solutions to the Yang- Baxter equation to obtain 
such representations. The section ends with a discussion of the role of the 
algebraic Yang-Baxter equation in configuring quantum circuit diagrams. 

In section 4, we work out details of the invariant of knots and links that 
is associated with the universal gate and give a number of examples. In 
particular, we show that this invariant measures the linking of the Borromean 
Rings and the Whitehead Link, both examples of links with zero linking 
numbers. 
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In section 5, we indicate how to formulate a quantum computation of a 
quantum link invariant in terms of a preparation, a unitary evolution and a 
measurement. We include in this context a process that quantum computes 
the absolute value of the trace of an arbitrary unitary transformation. These 
ideas are applied in Section 7 first to a unitary representation of the three 
strand braid group that will produce a good chunk of the Jones polynomial 
for three strand braids when configured as a quantum computer, and then 
to the invariant discussed in section 4. 

Entanglement is an integral ingredient in certain communications proce- 
dures such as the teleportation of quantum states. In section 6 we digress 
on the structure of teleportation, using the ideas prsented in the previous 
section for obtaining the trace of a unitary transformation. By associating 
a matrix M to a measurement state {M\ and using the entangled state \5) 
used for preparation in the trace calculation of the previous section, we show 
that for unitary M there is a full teleportation procedure for obtaining M\%1)) 
from a given state This discussion will be expanded in subsequent pa- 
pers to deal with the question of quantum computation in general, and the 
specific problem of computing knot invariants that are based on non-unitary 
solutions to the Yang-Baxter equation. The approach to teleportation given 
here is inherently topological (in the diagrammatic sense) and we shall take 
up its applications in subsequent papers [7]. 

In section 8 we discuss the relationship between topological entanglement 
and quantum entanglement. Wc recall an invariant of links associated with 
the solution to the Yang-Baxter equation used for Theorem 1. This solution, 
R' , makes an invariant that detects linking numbers of two-component links 
exactly when R' is capable of entangling quantum states. Examples like this, 
and invariants like the one constructed via the matrix i?, indicate relation- 
ships between topological entanglement and quantum entanglement. Other 
examples, such as the braid group representation representing the Jones poly- 
nomial of Section 4, do not exhibit such behaviour. The question remains 
open. In this section we give an example that effectively destroys the hope 
of continuing an analogy of Aravind that would identify the cutting of a link 
component with an observation of a state. Aravind himself showed that his 
notion was not invariant under basis change. We point out that it is easy to 
build states whose entanglement or lack of it after an observation is a mat- 
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ter of probability obtained from a probability amplitude. Since linking of 
classical links is not a matter of probability, this destroys the possibility of a 
direct relationship between classical linking and quantum entanglement. Of 
course, there may be more subtle avenues. We are in the process of working 
on such ideas. 

Section 9 is a capsule summary of topological quantum field theory from 
the point of view of anyonic models for quantum computation. We have 

included this section to indicate how braiding gates fit into a wider context. 
In section 10, we carry on a philosophical discussion about the relationship 
of quantum and topological entanglement, speculating that a spin network 
pregeometry of the right kind could enlighten us in our quest. 
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vanced Research Projects Agency (DARPA) and Air Force Research Labo- 
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and to give thanks to the University of Waterloo and the Perimeter Insti- 
tute in Waterloo, Canada for their hospitality during the preparation of this 
research. It gives both authors pleasure to thank Michael Nielson for very 
useful comments on an early version of this paper. We thank Ferando Souza, 
Hilary Carteret and Niel de Beaudrap for helpful conversations. 



2 Braiding Operators and Universal Gates 

We shall assume that the reader is familiar with the notions of knots, links 
and braids in Euclidean three dimensional space. Recall that a knot is an 
embedding of a circle, taken up to topological equivalence, and that a link is 
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an embedding of a collection of circles, taken up to topological equivalence. 
Braids form a group under concatenation, where the concatenation of two 
braids is obtained by attaching the bottom strands of the first braid to the 
top strands of the second braid. 

A class of invariants of knots and links called quantum invariants can 
be constructed by using representations of the Artin braid group, and more 
specifically by using solutions to the Yang-Baxter Equation [3] , first discov- 
ered in relation to 1 + 1 dimensional quantum field theory, and 2 dimensional 
statistical mechanics. Braiding operators feature in constructing representa- 
tions of the Artin Braid Group, and in the construction of these invariants 
of knots and links. 

A key concept in the construction of quantum link invariants is the as- 
sociation of a Yang-Baxter operator R to each elementary crossing in a link 
diagram. The operator i? is a linear mapping 

R: V®V — >V®V 

defined on the 2-fold tensor product of a vector space V", generalizing the 
permutation of the factors (i.e., generalizing a swap gate when V represents 
one qubit). Such transformations are not necessarily unitary in topological 
applications. It is a motivation for our research to understand when they can 
be replaced by unitary transformations for the purpose of quantum comput- 
ing. Such unitary i?-matrices can be used to make unitary representations 
of the Artin Braid group. 

A solution to the Yang-Baxter equation, as described in the last para- 
graph is a matrix i?, regarded as a mapping of a two-fold tensor product of 
a vector space V ®V to itself that satisfies the equation 

{R /)(/ R){R ®I)^{I®R){R® R). 

Prom the point of view of topology, the matrix it! is regarded as representing 
an elementary bit of braiding represented by one string crossing over another. 
In Figure 1 below, we have illustrated the braiding identity that corresponds 
to the Yang-Baxter equation. Each braiding picture with its three input lines 
(below) and output lines (above) corresponds to a mapping of the three fold 
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tensor product of the vector space V to itself, as required by the algebraic 
equation quoted above. The pattern of placement of the crossings in the 
diagram corresponds to the factors I and I <Si R. This crucial topological 
move has an algebraic expression in terms of such a matrix R. Our main ap- 
proach to relate topology, quantum computing, and quantum entanglement 
is through the use of the Yang-Baxter equation. In order to accomplish this 
aim, we need to study solutions of the Yang-Baxter equation that are unitary. 
Then the it! matrix can be seen either as a braiding matrix or as a quantum 
gate in a quantum computer. 




Figure 1 The Yang-Baxter Equation - 

{R®I){I®R){R®I) = {I®R){R®I){I®R). 
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Braid Generators 



s^-1s^ =1 



^1 S2 s-j = S2 s-j S2 



3 ^1 



Figure 2 - Braid Generators and Relations 

The problem of finding solutions to the Yang-Baxter equation that are 
unitary turns out to be surprisingly difficult. Dye [12] has classified all such 

matrices of size 4 x 4. A rough summary of her classification is that all 
4x4 unitary solutions to the Yang-Baxter equation are similar to one of the 
following types of matrix: 



I l/\/2 l/^/2 \ 

1/^2 -1/^2 

1/^ l/^h 

-l/\/2 1/^2 I 



9 



R' = 
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where a,b,c,d are unit complex numbers. 

For the purpose of quantum computing, one should regard each matrix 
as acting on the stamdard basis {|00), |01), |10), |11)} of H = V ^V, where 
y is a two-dimensional complex vector space. Then, for example we have 

it:|oo) = (i/v^)|oo)-(i/^)|ii), 

i?|01) = (l/v^)|01) + (l/v^)|10), 

i?|10) = -(1/V2)|01) + (1/V2)|10), 

= (l/v^)|00) + (1/72)111). 

The reader should note that R is the familiar change-of-basis matrix from 
the standard basis to the Bell basis of entangled states. 

In the case of i?', we have 

R'\00) = a|00),i?'|01) = c|10), 

R'\10) ^b\01),R'\ll) = ci|ll). 

Note that R' can be regarded as a diagonal phase gate P, composed with a 
swap gate S. 

/ a \ 

6 
c 
V d / 
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s = 



( 1 








M 








1 








1 








I 








1 / 



Compositions of solutions of the (Braiding) Yang-Baxter Equation with the 
swap gate S are called solutions to the algebraic Yang-Baxter equation. Thus 
the diagonal matrix P is a solution to the algebraic Yang-Baxter equation. 



2.1 Universal Gates 

A two-qubit gate G is a unitary linear mapping G : V (E)V — > V where V is 
a two complex dimensional vector space. We say that the gate G is universal 
for quantum computation (or just universal) if G together with local unitary 
transformations (unitary transformations from V to V) generates all unitary 
transformations of the complex vector space of dimension 2" to itself. It is 
well-known [44] that CNOT is a universal gate. 

A gate G, as above, is said to be entangling if there is a vector 

\a(3) ^\a)®\(3) eV®V 

such that G\aj3) is not decomposable as a tensor product of two qubits. 
Under these circumstances, one says that G\a(3) is entangled. 

In [6], the Brylinskis give a general criterion of G to be universal. They prove 
that a two-qubit gate G is universal if and only if it is entangling. 

The reader will also be interested in the paper [5] and the url 

http : / /www. physics. uq.edu. au/ gqc/ , 

wherein the practical algorithm in [5] , for expressing entangling gates in terms 
of CNOT and local transformations, is implemented online. 
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It follows at once from the Brylinski Theorem that the matrices R, R', and 
R" are universal gates, except for certain specific choices of parameters in 
R' and R". In a sequel to this paper [13] we will give a complete catalogue 
of universality for two-qubit gates that are solutions to the Yang-Baxter 
equation. In this paper, we shall concentrate on specific examples and their 
properties. 

Remcirk. A two-qubit pure state 

10) = a|00) + 6|01) + c|10) + d\ll) 

is entangled exactly when {ad — be) ^ 0. It is easy to use this fact to check 
when a specific matrix is, or is not, entangling. 

Theorem 0. Let D denote the phase gate shown below. D is a solution to 
the algebraic Yang-Baxter equation (see the earlier discussion in this section). 
Then D is a universal gate. 



D = 



/ 1 













1 














1 





V 








-1 / 



Proof. It follows at once from the Brylinski Theorem that D is universal. 
For a more specific proof, note that CNOT = QDQ^^, where Q — H ® I , 
H is the 2x2 Hadamard matrix. The conclusion then follows at once from 
this identity and the discussion above. We illustrate the matrices involved 
in this proof below: 



« = (1/^/2) ( ; _\ ) 



Q = (1/^2) 



/ 1 


1 
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-1 
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I 





1 
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/ 1 








\ 





1 














1 





V 








-1 ) 



QDQ-^ = QDQ 



This completes the proof of the Theorem. 



/ 1 













1 
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V 





1 


/ 



CNOT 



□ 



Remcirk. We thank Martin Roetteles [46] for pointing out the specific fac- 
torization of CNOT used in this proof. 

Theoreml. The matrix solutions R' and R" to the Yang-Baxter equation, 
described above, are universal gates exactly when ad— be ^ for their internal 
parameters a, b, c, d. In particular, let /?o denote the solution R' (above) to 
the Yang-Baxter equation with a = b = c = l.d = —1. 

/ 1 \ 
10 
1 




Rq 





V 










-1 / 



Then Rq is a universal gate. 



Proof. The first part follows at once from the Brylinski Theorem. In fact, 
letting H be the Hadamard matrix as before, and 



1/V2 i/V2 \ X / 1/^2 1/V2 \ 



{l-t)/2 {l + t)/2 
(1-0/2 (-l-i)/2 



Then 



CNOT = (A (8) pi){Ro{I <8) cr)Ro){H H). 

This gives an explicit expression for CNOT in terms of Rq and local unitary 
transformations (for which we thank Ben Reichardt in response to an early 
version of the present paper). □ 
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Remark. Let SWAP denote the Yang-Baxter Solution R' with a 
d=l. 



SWAP 



SWAP is the standard swap gate. Note that SWAP is not a universal gate. 
This also follows from the Brylinski Theorem, since SWAP is not entangling. 
Note also that Rq is the composition of the phase gate D with this swap gate. 



/ 1 








\ 








1 








1 








I 








1 J 



Theorem2. Let 



R 



( 1/V2 







1/^2 -l/x/2 

1/^2 1/v^ 





1/V2 \ 





be the unitary solution to the Yang-Baxter equation discussed above. Then 
is a universal gate. The proof below gives a specific expression for CNOT 
in terms of R. 

Proof. This result follows at once from the Brylinksi Theorem, since R is 
highly entangling. For a direct computational proof, it suffices to show that 
CNOT can be generated from R and local unitary transformations. Let 

' 1/v^ \ 
^ 1/^2 -1/72 ] 



a 



7 



-l/^/2 1/V2 

i/V2 i/V2 

1/V2 i/V^ 
1/V2 -i/V2 



-1 

Let M — a<Si P and N — ^ ^ S. Then it is straightforward to verify that 

CNOT = MRN. 

This completes the proof. □ 
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Remark. Wc take both Theorems 1 and 2 as suggestive of fruitful inter- 
actions between quantum topology and quantum computing. It is worth 
comparing these Theorems with the results in [16], a comparison that we 
shall leave to a future paper. 

Remark. We thank Stephen Bullock for his help in obtaining this result. 
On showing him the Yang-Baxter solution R used in the above proof, he 
showed us the paper ([48] by V. V. Shende, S. S. Bullock and I. L. Markov) 
in which he and his coauthors give a criterion for determining if a 4 x 4 
unitary matrix can be generated by local unitary transformations and a sin- 
gle CNOT. We then calculated that criterion and found that R passes the 
test. Bullock then showed us how to apply their theory to obtain the specific 
transformations needed in this case. Thus the above result is a direct ap- 
plication of their paper. The criterion also shows that the solutions of type 
R and R listed above require two applications of CNOT. We will discuss 
their structure elsewhere; but for the record, it is of interest here to record 
the Shende,Bullock, Markov criterion. 

Theorem [48]. We shall say that a matrix can be simulated using k CNOT 
gates if it can be expressed by that number k of CNOT gates plus local 
unitary transformations. Let E be the following matrix: 



E^ 



1 ° 








M 








-1 








-1 








1 1 








I 



Let [/ be a matrix in SU{A). Let ^{U) be defined by the formula 

7(t/) = UEUE. 
Let tr{M) denote the trace of a square matrix M. Then 

1. U can be simulated using zero CNOTS if and only if 'y{U) = I, where 
/ denotes the identity matrix. 

2. U can be simulated using one CNOT gate if and only if tr[j{U)] — 
and ^{Uy = -I. 
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3. U can be simutated using two CNOT gates if and only if tr[7(t/)] is 
real. 

□ 

Note that in applying this criterion, the matrix in question must be in the 
special unitary group. We leave it to the reader to show that matrices of 
type R' and R" require two CNOTS, and that the matrix R is picked by 
this criterion to require only one CNOT, just as we have shown explicitly 
above. Note that since R^ is the identity, we have R~^ = R^, as well as the 
fact that R^^ can be expressed in terms of local transformations and a single 
application of CNOT. 

3 Generalizing and Representing the Artin 
Braid Group 

Let Bn denote the Artin braid group on n strands [45]. We recall here that 
Bn is generated by elementary braids {si, • • • , with relations 

1. SiSj = SjSi for |i — j| > 1, 

2. SiSi+iSi = Si+iSiSi+i for i = 1, • • - n - 2. 

See Figure 2 for an illustration of the elementary braids and their relations. 
Note that the braid group has a diagrammatic topological interpretation, 

where a braid is an intertwining of strands that lead from one set of n points 
to another set of n points. The braid generators Sj are represented by dia- 
grams where the i-th and {i + l)-th strands wind around one another by a 
single half-twist (the sense of this turn is shown in Figure 2) and all other 
strands drop straight to the bottom. Braids are diagrammed vertically as 
in Figure 2, and the products are taken in order from top to bottom. The 
product of two braid diagrams is accomplished by adjoining the top strands 
of one braid to the bottom strands of the other braid. 
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In Figure 2 we have restricted the illustration to the four-stranded braid 
group B4. In that figure the three braid generators of B4 are shown, and 
then the inverse of the first generator is drawn. Following this, one sees the 
identities SiSi^ — 1 (where the identity element in B4 consists in four vertical 
strands), S1S2S1 = S2S1S2, and finally S1S3 = S3S1. With this interpretation, 
it is apparent from Figures 1 and 2 that the second braiding relation (above) 
is formally the same as the Yang-Baxter equation. 

In fact, if V denotes the basic vector space (the space for one qubit in 
our context), and R is an invertible solution to the Yang-Baxter equation as 
described in section 2, then we obtain a representation of the n strand braid 
group into the vector space of automorphisms of the n-th tensor power of V : 

repn : A«t(y^") 

by defining 

repnisi) = /®'-^ (g) 

where I denotes the identity mapping on V. Note that since i? is a unitary 
matrix, it follows that this representation is unitary. We shall call this the 
standard method for making a representation of the braid group from an 
invertible solution to the Yang-Baxter equation. Note that repn{si) is sup- 
ported by R on the i-th and {i + l)-th tensor factors and is the identity 
mapping on the other factors. 

We now wish to generalize the classical Artin braid group to a larger group 
whose representations can include compositions of the elements repn{si) con- 
structed in the last paragraph with local unitary transfomations. A diagram- 
matic example of such a composition is given in Figure 3, where we illustrate 
(as in Theorem 2) the expression of CNOT in terms of R and local unitary 
transfomations. In this diagram, the local unitary transformations are ini- 
cated by nodes on single braiding lines. This corresponds to the fact that the 
local unitary operations act on single tensor factors. Thus we shall generalize 
the braid group so that the new group generators arc represented on single 
tensor factors while the elementary braids are represented on two essential 
tensor factors. This is formalized in the next paragraph. 

Let G be any group, and let G"^" denote the n-fold tensor product of G 
with itself, where by this tensor product we mean the group whose elements 



17 



are of the form (71 ® (72 C?) ■ ■ ■ ® (?„. with gi G G, satisfying 

(51 (g) 5(2 (8) • • • (g) 5f„)(5fi (g) 52 ® • • • O g'n) = 5151' <8) 9292 • • • S-nSn'- 

We articulate G®" in tensor language, because we wish to consider represen- 
tations of the group where individual members of G go to matrices and the 
elements 91® 92® ■ ■ ■ ® 9n are sent to tensor products of these matrices. Let 
hi{9) = e®e---®9®e---®e, where e is the identity element in G, and the 
element 9 is in the i-th place in this tensor product string. Then 

91® 92® ■ ■ ■ ® 9n = hi{gi) ■ ■ ■ hn{gn), 

and G®" is generated by the elements hi{g) where i ranges from 1 to n and 
g ranges over elements of G. Note that for i j, and for any g,g' e G, the 
elements hi{g) and hj{g') commmute. Note that hi is an isomorphism of G 
to the subgroup hi{G) of G®". 

We define an extension GBn of the braid group S„ by the group G®" as 
follows: GBn is freely generated by G®" and En modulo the relations 

hi{9)sj = Sjhi{g) 

for all g' in G and all choices of i and j such that i<j or i>j + l. 

Just as there is a diagrammatic intepretation of the braid group Bn in 
terms of strings that entangle one another, there is a diagrammatic intepre- 
tation of GBn- Think of a braid diagram and suppose that on the lines of 
the diagram there is a collection of labelled dots, with each dot labeled by an 
element of the group G. Let it be given that if two dots occur consecutively 
on one of the strings of the braid diagram, then they can be replaced by a 
dot labeled by the product of these two elements. We make no assumptions 
about moving dots past elementary braiding elements (which have the ap- 
pearance of one string passing over or under the other). It is easy to see 
that this diagrammatic description also defines a group extending the braid 
group, and that this diagrammatic group is isomorphic to GBn- 

We apply this description of G-B„ by taking G = U{2), the 2x2 unitary 
matrices, viewing them as local unitary transformations for quantum com- 
puting. Then we let UBn denote U{2)Bn and take the representation of UBn 
to Aut(y®'^) that is obtained by the mapping 

r : UBn — > AutiV®'') 
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defined by T{hi{g)) = hi{g) for g in U{2) and T{si) = I^'-^ ® i? ® jm-i-i 
where R is the Yang-Baxter solution discussed in Theorem 2. 

Recall that a quantum computer is an n-qubit unitary transformation 
U coupled with rules/apparatus for preparation and measurement of quan- 
tum states to which this transformation is applied. We then conclude from 
Theorem 2 that 



Theorem 3. Any quantum computer has its basic unitary transformation 
U in the image of F. 




R CNOT 
CNOT =(A(g)B)R(C(g)D) 

Figure 3 - CNOT ^(A® B)R{C ® D) 

Remark. This theorem means that, in principle, one can draw a circuit dia- 
gram for a quantum computer that is written in the language of the extended 
braid group UBn. In particular this means that braiding relations will apply 
for sectors of the circuitry that are not encumbered by local unitary transfor- 
mations. Typically, there will be many such local unitary transformations. 
We will investigate this braid algebraic structure of quantum computers in 
a sequel to this paper. A key illustration of Theorem 3 is the diagrammatic 
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interpretation of Theorem 2. This is shown in Figure 3 where we have written 
in diagrams an equation of the form 

CNOT ^{A® B)R{C (g) D). 

Here A, C and D represent the local unitary matrices that are used in 
Theorem 2 (with different names) to express CNOT in terms of R. 

In general, a unitary transformation can be written as an extended braid- 
ing diagram, with appearances of R occupying two adjacent strands, and of 
local transformations occupying single strands. Note that since R^ is the 
identity. Theorem 3 actually says that a unitary transformation can be built 
via a representation of the extension of a quotient B'^ of the braid group S„, 
where each braid generator s- in B'^ has order eight. It is worth investigating 
the algebraic structure of B'^. This is a topic for further research. 



3.1 The Algebraic Yang-Baxter Equation 

If R denotes a solution to the Yang-Baxter equation (not neccessarily the 
R of Theorem 2), then we can consider the composition r = SR where S 
is the swap gate defined in section 2. If we think of r as supported on two 
tensor lines, and write rij for the same matrix, now supported on tensor lines 
i and j, (all other lines carrying the identity matrix) then we find that the 
Yang-Baxter equation for R is equivalent to the following equations for Vij. 

The above equation is called the Algebraic Yang-Baxter Equation. See Figure 

4 for an illustration of this relationship. In making circuit diagrams to apply 
Theorem 1, it is useful to use the formalism of the algebraic Yang-Baxter 
equation since we can then think of the phase gate D as such a solution and 
use the above relation to relocate compositions of D on diffferent tensor hues. 
Given a solution to the algebraic Yang-Baxter equation, plus the swap gate 

5 we can again define a generalization of the Braid Group that includes local 
unitary transformations on the single tensor lines. We will leave detailed 
application of this point of view to a sequel to this paper. 
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Figure 4 -Algebraic Yang Beixter Equation 

4 An Invariant of Knots and Links Associ- 
ated with the Matrix R 

A well-known relationship between braids and knots and links allows the 
construction of invariants of knots and links from representations of the Artin 
Braid Group. We give here a quick summary of these relationships and then 
apply them to the quantum universal matrix i?, showing that it gives rise to 
an interesting invariant of knots and links. The reader should note that this 
section is concerned only with the classical braid group. It does not use the 
extensions of the braid group that are discussed in the previous section. 
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At this point it is worth making a digression about the Reidemeister 
moves. In the 1920 's Kurt Reidemeister proved an elementary and important 
theorem that translated the problem of determining the topological type of 
a knot or link to a problem in combinatorics. Reidemeister observed that 
any knot or hnk could be represented by a diagram, where a diagram is a 
graph in the plane with four edges locally incident to each node, and with 
extra structure at each node that indicates an over-crossing of one local arc 
(consisting in two local edges in the graph) with another. See Figure 5. The 
diagram of a classical knot or link has the appearance of a sketch of the 
knot; but it is a rigorous and exact notation that represents the topological 
type of the knot. Reidemeister showed that two diagrams represent the same 
topological type (of knottedness or linkedness) if and only if one diagram 
can be obtained from another by planar homeomorphisms coupled with a 
finite sequence of the Reidemeister moves illustrated in Figure 6. Each of 
the Reidemeister moves is a local change in the diagram that is applied as 
shown in this Figure. 




Figure 5 - A Knot Diagram 
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Figure 6 - Reidemeister Moves 




We say that two knots or links are isotopic if one can be obtained from 
the other by a sequence using any of the three Reidemeister moves (plus 
global topological mappings of the diagram plane to itself). 

The first significant fact relating hnks and braids is the 

Theorem of Alexander. Every knot or link is isotopic to the closure of a 
braid. 

The closure of a braid b, here denoted CL{b), is obtained by attaching each 
top strand to the corresponding bottom strand in the fashion shown in Figure 
7. The closed braid is a weave that proceeds circularly around a given axis. 
There are many proofs of Alexander's Theorem. The interested reader should 
consult [45]. 




Hopf Link 




Trefoil Knot 




igure Eight Knot 



Figure 8 - Closing Braids to Produce Hopf Link, Trefoil Knot and 

Figure Eight Knot 
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Given that every knot or link can be represented by a closed braid, it is 
natural to wonder whether the classification of braids will effect a classifica- 
tion of the topological types of all knots and links. The situation is more 
complicated than one might have expected. There are many braids whose 
closure is isotopic to any given knot or link. Here are two basic methods for 
modifying a braid b in so that the topological type of its closure does not 
change: 

1. Let g be any braid in B^- Then 

CL{gbg-') = CL{b) 

where we use equality to denote isotopy of knots and links as described 
above. 

2. Note that if b is in then bsn is in Bn+i- It is easy to sec that 

CL{b) = CL{bSn) 

and 

CL{b)=CL{bs-'). 

In light of the equvalences we have just indicated, the following two moves 
on braids are called the Markov moves (after Markov who enunciated the 
Theorem we state below): 

1. Markov Move 1 Replace a braid b by gbg~^ where g is another braid 
with the same number of strands. 

2. Markov Move 2 Replace a braid b E B^hy either 6s„ or by bs~^ or 
vice versa, replace bs^^ with b. 
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CL( b s 3 ) = CL(b) 



Figure 9 - Illustration of the Second Markov Move 
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tr(rep4(bS3)) = sqrt{2) Mrepg (b)) 



Figure 10 - Illustration of the Behaviour of the Trace on the 

Second Markov Move 
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Markov Theorem. Suppose that b and b' are two braids (of possibly differ- 
ent numbers of strands) with CL{b) = CL{b'). Then b' can be obtained from 
6 by a series of braid equivalences coupled with applications of the Markov 
moves) . 

Remcirk. For proofs of the Markov Theorem, see [4] and [33]. See Figure 
9 for an illustration of the second Markov move. Notice that in making this 
move wc promote the braid 6 G i?n to a braid in by adding a right-most 
strand. Then we multiply by s„ G -Bn+i. The closure of the resulting braid 
differs by a single first Reidemeister move from the closure of b. The upshot 
of this theorem is that it is possible for a trace function on a representation 
of the braid group to give rise to topological information about the closure 
of the braid. For example, suppose that we have a unitary representation of 
the braid group arising from a unitary solution of the Yang-Baxter equation, 
as described in Section 3. Let the representation be denoted by 

Let 

T{b) = tr{repn{b)) 

where tr{M) denotes the trace of a square matrix M. Then, since the trace 
of any linear mapping satisfies tr{AB) = tr{BA), it follows that T{gbg~^) — 
t(6), and hence r gives the same values on braids that differ by Markov 
moves of type 1. We would like r to be invariant under Markov moves of 
type 2, but this is usually too much to ask. It is standard practice in the 
literature of link invariants to search for a matrix rj mapping V to V such 
that the modified trace TR{b) = tr{r]'^"'repn{b)) has a multiplicative property 
under the second Markov move in the sense that TR{bSn) — aTR{b) and 
TR{bs~^) — a~^TR{b), where o; is a invertible constant in the ring of values 
for the trace. Such a function TR is called a Markov trace, and one can 
normalize it to obtain a function that is an invariant of isotopy of links 
by defining 1(b) = a~^'^^^TR{b), where w{b) is the sum of the signs of the 
crossings of the braid b. 

In the case of our computationally universal matrix /?, the bare trace 
r(6) = tr{repn{b)) behaves in a very simple way under the second Markov 
move. We find (and will show the details below) that 

T{bsn) = \/2r(6) 
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and 

T{bs-') = V2T{b). 

Note that the multiphcative factor is the same for both types of second 
Markov move. Instead of making a normahzing factor from this, we can say 
that if two hnks CL{b) and CL{b') are isotopic, then t(6) and t(6') will differ 
by a multiplicative factor that is some power of the square root of two. In 
particular, this means that if r(6) and r(6') have different sign, or if one is 
zero and the other not zero, then we know that the closures of b and b' are 
not isotopic. 

Here is the matrix R. 

/ l/\f2 1/^2 \ 

1/v^ -1/V2 

1/V2 1/V2 
^ -l/\/2 l/\/2 J 

We have 

R\ab) = R'^lm + Rll\Ql) + Rl^\lQ) + Rll\ll). 

Let 

tr^iR) = EkKt 

be the partial trace of R with respect to the second tensor factor. Then 

tr2{R)^{l/V2 + l/V2)I ^V2I 

where / denotes the 2x2 identity matrix and 

tr2{R-^)^V2I. 

Recall from the previous section that 

repn : 5„ — > Aut{V®'') 

is defined on braid generators by the equation 

repnisi) = I^'-^ i? ® 

where / denotes the identity mapping on V. Now suppose that b G Bn- We 
want to compute T{bsn). This expression requires interpretation. When we 
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write bsn, we are taking b in Bn and regarding it as an element of Bn+i by 
adding an extra right-most strand to b. In general, by adding strands in this 
way, we have standard embeddings of S„ in when m > n. Working in 
Bn+i, we have 

repn+i{b) = repn{b) (g) / 

and 
Thus 

T{bsn) = tr{repn+i{bsn)) = tr{{rep{b) ® ® R)). 

From this it is easy to see that in tracing rep„_|_i(6s„), the rightmost indices 
of the matrix R (in the n + 1 tensor factor) are contracted directly with one 
another (since b is supported on the first n strands). Thus the partial trace 
is applied to the R that appears in the representation of bsn corresponding 
to Sn- It follows from this that 

T{bsn) = tr{repn{b))tr2{R) = r(6)V2. 

Hence, for the tensor representation built from R as described in section 3, 
we have 

r{bsn) = V2T{b) 

and in like manner, we have 

T{bs-')^V2T{b). 

This proves the assertions we have made about the properties of r for this 
R. 

Remcirk. In Figure 10, we illustrate diagrammatically the above argument 
at the index level. In this illustration, we have placed a shaded box around 
a braid to indicate the application of the representation of the braid group. 
Thus a shaded braided box represents a matrix with upper indices corrcp- 
sponding to the upper strands on the box, and lower indices corresponding 
to the lower strands on the box. The simplest instance of such a matrix is a 
single vertical line which represents the identity matrix, and iconically indi- 
cates the identity of the top index with the bottom index (hence representing 
the identity as a Kronecker delta). 
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It is a fact that shaded boxes so placed on the braids give a correct picture 
of the contractions of the corresponding matrices via the convention that we 
contract the indices along lines that connect free index ends between diagram- 
matic matrices. The figure then illustrates directly via these diagrammatic 
matrices how we obtain the formula 

tr{repn+iibsn)) = \/2tr{repn{h)). 

Note that the trace of a diagrammatic matrix has exactly the same form as 
the closure of a braid, since the connection of two open lines corresponds to 
the identification and contraction over their respective indices. 

Finally, here is the same argument using matrix algebra with indices. We 
use the Einstein summation convention: Summation is taken over repeated 
upper and lower indices. Suppose that repn{b) = (M^^*) where a and b are 
vectors of indices for the the first n — 1 factors of the tensor product, and i 
and j are individual indices with values or 1. Then 

where is the 2x2 identity matrix. Furthermore 

repn+i{sn) = S^Rl;%. 

Hence 

from which it follows that 

trirepn^bsn)) = tr{M^:;i^ = tr{M^:;^Q 
where in the last equality we have contracted the last indices of R. Since 

it follows that 

tr{repn+i{bsn)) = tr{M'^;;V28i) - V2tr{M^;l) = V2tr{repn{b)). 

This completes the explicit index verification of the behaviour of r under the 
second Markov move. 
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This invariant of knots and links turns out to be quite interesting. For 
example, it detects the linkedness of the Borromean rings (depicted in Figures 
7 and 16). It gives the following values: 

1. riUnlink of Three Components) = 8 > 

2. T{Hopf Link) = 

3. r{Trefoil Knot) = -2^2 < 

4. T{Figure Eight Knot) = -4 < 

5. T{Borromean Rings) = — 8 < 

Note that r does not detect the difference between the trefoil knot, the 
figure eight knot and the Borromean rings, but it does show that the Hopf 
link is linked, that the Borromean rings are linked, and that the trefoil knot 
and the figure eight knot arc knotted. See Figure 8 for illustrations of these 
knots and links. It remains to be seen how the quantum entangling properties 
of the matrix R are related to the behaviour of this link invariant. 

Remark on a Skein Relation. In this subsection, we point out that there 
is a skein relation that helps in the computation of the trace r(6) for a braid h. 
A skein relation is an equation about an invariant involving local changes at 
the site of a single crossing in corresponding braid or link diagrams. The first 
skein relation in knot theory was discovered and utilized for the Alexander 
polynomial by John H. Conway in his remarkable paper [8]. (Conway used 
an idea that was implicit in Alexander's original paper of 1928.) The Jones 
polynomial and many other knot polynomial invariants satisfy such relations. 

The matrix R satisfies the equation 

R + R-' = V2I2 

where /„ denotes the 2" x 2" identity matrix. We leave it to the reader to 
check this fact. It is also easy to check that 

R' = l2 
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and that all the lower powers are non-trivial. The fact that R has finite order 
certainly limits its power as a link or braid invariant. For example, we have 
the eight-fold periodicity 

as a direct consequence of the finite order of R. On the other hand the identity 
R + R-^ = V2I2 can be viewed as a method for simplifying the calculations 
for a braid. This implies the skein relation 

r(b)+r{b') = V2r(b") 

when b and b' are elements of the n-strand braid group that differ at a single 
crossing and b" is the result of replacing this crossing by an identity braid. 
The crossing can be interpreted as a single instance of Sj for some i, and we 
then use repn{otSi (3) + repn{as~^ P) — y/2 repn{oiP). 

Example 1. Here is the simplest example of this sort of computation. We 
work in B2 and let s = Si : 

Here we have r{ss-^) = r(/2) = 4 and r(s) = V2r(/i) = 2^2. Hence 
t{ss) = 0, as we remarked earlier with ss the braid representative for the 
Hopf link (See Figure 8). More generally, we have 

r(s"+i) +r(s"-^) = V2r(s") 

so that 

Letting 1 denote the identity braid in B2, we then have 

r(l)=4 

t(s) = 2V2 
t(s') = 
r(s^) = -2V2 
t(s^) = -4 
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-2\/2 



r(s^) = 
r(/) = 2V2 

with the periodicity 

r(s"+^) = r(s"). 

Since s'^ and s"^ close to knots, we see that this invariant can distinguish these 
two knots from one another, but cannot tell that the closure of is knotted. 

Example 2. Let 6 = sls2^ siS2^ ■ See Figure 11. The closure of b is W, a 
link of two components, with linking number equal to zero. W is called the 
Whitehead Link, after the topologist, J. H. C. Whitehead, who first studied 
its properties. We shall check that r(6) = — 4\/2, showing that our invariant 
detects the linkedness of the Whitehead link. 

We use skein relation for the first appearance of from the left on the 
word for b. This gives 

r(6) = -t{sIs2SiS2^) + V2t{sIsiS2^). 

Note that 

sls2SiS2^ = Sl{siS2Si)s2^ = Sl(s2SlS2)s2 ^ = S1S2S1. 

Then 

t(siS2Si) = t((siS2)si) = t(si(siS2)) = r{sls2) = V2t{sI) = 0. 

Hence 

r(6) = V2t{sIsiS2^) = V2r(s?s^^) = V2V(s?) = 2{-2V2) = -4^2. 
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Figure 11 - Whitehead Link, W = CL{b = sfs^^sis^^). 

5 Quantum Computation of Knot Invariants 

Can the invariants of knots and links such as the Jones polynomial be con- 
figured as quantum computers? This is an important question because the 
algorithms to compute the Jones polynomial are known to be A^P-hard, 
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and so corresponding quantum algorithms may shed hght on the relation- 
ship of this level of computational complexity with quantum computing (See 
[19]). Such models can be formulated in terms of the Yang-Baxter equation 
[22, 23, 28, 25]. The next paragraph explains how this comes about. 

In Figure 12, we indicate how topological braiding plus maxima (caps) 
and minima (cups) can be used to configure the diagram of a knot or link. 
This also can be translated into algebra by the association of a Yang-Baxter 
matrix R (not necessarily the R of the previous sections) to each crossing 
and other matrices to the maxima and minima. There are models of very 
effective invariants of knots and links such as the Jones polynomial that 
can be put into this form [25]. In this way of looking at things, the knot 
diagram can be viewed as a picture, with time as the vertical dimension, of 
particles arising from the vacuum, interacting (in a two-dimensional space) 
and finally annihilating one another. The invariant takes the form of an 
amphtude for this process that is computed through the association of the 
Yang-Baxter solution R as the scattering matrix at the crossings and the 
minima and maxima as creation and annihilation operators. Thus we can 
write the amplitude in the form 

Zk = {CUP\M\CAP) 

where {CUP\ denotes the composition of cups, M is the composition of ele- 
mentary braiding matrices, and \CAP) is the composition of caps. We regard 
{CUP\ as the preparation of this state, and \CAP) as the measurement of 
this state. In order to view Zk as a quantum computation, M must be a 
unitary operator. This is the case when the i?-matrices (the solutions to 
the Yang-Baxter equation used in the model) are unitary. Each i?-matrix is 
viewed quantum gate (or possibly a composition of quantum gates), 

and the vacuum-vacuum diagram for the knot is interpreted as a quantum 
computer. This quantum computer will probabilistically (via quantum am- 
plitudes) compute the values of the states in the state sum for Zk- 
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{CAP\ [M easurement) 




Quantum Computation 



Unitary Braiding 



\CUP) (Preparation) 



Zk = {CAP\M\CUP) 



Figure 12 A Knot Quantum Computer 

The form of the model proposed for translating the Jones polynomial 
to a quantum computation is also the form of models for anyonic quantum 
computation (See [14, 15, 16, 17, 18]). In an anyonic model, the braiding 
corresponds to the motion of configurations of particles in a two-dimensional 
space. These theories are directly related to quantum link invariants and 
to topological quantum field theories [50]. It is hoped that quantum com- 
puting placed in the anyonic context can be made resistant to the effects of 
dccohcrence due, in part, to the invariance of topological structures under 
perturbation. 

The formalism of configuring a computation in terms of preparation and 

measurement in the pattern of Figure 12 can be used in very general quantum 
computational contexts. For example, let U he a, unitary transformation 
on H = V^"^ where V is the complex two-dimesional space for a single 
qubit. Represent U as a box with n input lines at the bottom and n output 
lines at the top, each line corresponding to a single qubit in an element of 
the tensor product H with basis {|q;)|q; is a binary string of length n}. Let 
\S) = Sct|«,a) G H ® H where a runs over all binary strings of length 
n. Note that (5| is the following covector mapping H ® H to the complex 
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numbers C : 

{S\a, — 1 if a — P and {S\a, (3) — otherwise. 

Now let W — U <Si Ih: where In denotes the identity transformation of H to 
H. Then 

{6\W\6) = {6\U^Ih\6) = 

{6\^^U2\l,a) = i:^U:=tr{U). 

For example, {S\S) — 2^ — trilu). See Figure 13 for an illustration of this 
process. 



r 



(5| 



V ®1 



H 



J 



\5) 



Figure 13 - A quantum process to obtain \tr{U)\. 



Thus we see that we can, for any unitary matrix [/, produce a quantum 
computational process with preparation \5) and measurement {5\ such that 
the amplitude of this process is the trace of the matrix U divided by (\/2)". 
This means that the corresponding quantum computer computes the proba- 
bility associated with this amplitude. This probability is the absolute square 
of the amplitude and so the quantum computer will have |tr(C/)p/2" as the 
probability of success and hence one can find \tr{U) \ by successive trials. We 
have proved the 
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Lemma. With the above notation, the absolute value of the trace of a 
unitary matrix U, \tr{U)\, can be obtained to any desired degree of accuracy 
from the quantum computer corresponding to U <S> Ir with preparation \S) 
and measurement {S\, where \S) = Eq,|q;, a) e H ® H. 

The proof of the Lemma is in the discussion above its statement. 

6 Unitary Representations and Teleportation 

The formalism we used at the end of the last section to describe the (absolute 
value) of the trace of a unitary matrix contains a hidden teleportation. It is 
the purpose of this section to bring forth that hidden connection. 

First consider the state 

\5) = T.oc\oi,a) e H ^H. 

from the last section, where H = V^"' and K is a single-qubit space. One 
can regard \5) as a generahzation of the EPR state (00 1 + (11 1. 

Let G -ff be an arbitrary pure state in H. Let be an abitrary 
element of the dual of H H and consider the possibility of a successful 
measurement via in the first two tensor factors of 

The resulting state from this measurement will be 

{M\m5)]- 

If 

(A^| = E,,^M,,^(a|(/?|, 

then 

{M\[m5)] = S,,/5M,,/3(«|(/3|S,>,|7)|A)|A) 
= S«,^M,,^E^>^(a|7)(/3|A)|A) 

Thus we have proved the 
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Teleportation Lemma. Successful measurement via (A^| in the first two 

tensor factors of 



results in the state M'^\ip) where the matrix M represents the measurment 
state {M.\ in the sense that 



and denotes the transpose of the matrix M. 

This Lemma contains the key to teleportation. Let |^) be a state held by 
Ahce, where Alice and Bob share the generahzed EPR state \S). Ahce mea- 
sures the combined state reports to Bob that she has succeeded in 
measuring via (Ai \ (from some list of shared transformations that they have 
in common) by a classical transmission of information. By the Lemma, Bob 
knows that he now has access to the state M'^\ip). In this generalized version 
of teleportation, we imagine that Alice and Bob have a shared collection of 
matrices M, each coded by a bit-string that can be transmitted in a classical 
channel. By convention, Alice and Bob might take the zero bit-string to de- 
note lack of success in measuring in one of the desired matrices. Then Alice 
can send Bob by the classical channel the information of success in one of 
the matrices, or failure. For success. Bob knows the identity of the resulting 
state without measuring it. See Figure 13.1 for a schematic of this process. 




M|(p> 



<M| 
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Figure 13.1 - Matrix Teleportation 



In the case of success, and if the matrix M is unitary, Bob can apply 
{M'^)~^ to the transmitted state and know that he now has the original 
state lip) itself. The usual teleportation scenario, is actually based on a 
list of unitary transformations sufficent to form a basis for the measurement 
states. Lets recall how this comes about. 



First take the case where M is a unitary 2x2 matrix and let (71,(72,(73 
be the three Pauli matrices 





" 1 " 




" 


1 " 




" i' 




-1 


,Cr2 = 


1 





,0-2 = 


-i 



We replace (73 by —ia^ (for ease of calculation) and obtain the three matrices 
X,Y,Z: 



X 



1 
-1 



Y 



1 

1 



1 
-1 



Basis Lemma. Let M be a 2 x 2 matrix with complex entries. Let the 
measuring state for M be the state 

{M\ = MoolOO) + MoilOl) + MiollO) + Mn|ll). 

Let {XAi I denote the measuring state for the matrix XM (similarly for YM 
and ZM). Then the vectors 

{{M\,{XM\,{yM\,{ZM\} 

are orthogonal in the complex vector space V ^ V ii and only if M is a 
multiple of a unitary matrix U of the form 



U 



z w 
—w z 



with complex numbers z and w as generating entries. 

Proof. We leave the proof of this Lemma to the reader. It is a straightfor- 
ward calculation. 
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This Lemma contains standard teleportation procedure when one takes 
M = / to be the identity matrix. Then the four measurement states 

m{x\,{y\,{z\} 

form an orthogonal basis and by the Telportation Lemma, they successfully 
transmit {\'^),X'^\'^),Y'^\'^),Z'^\'^)} respectively. Bob can rotate each of 
these received states back to by a unitary transformation (Remember 
that states are determined up to phase.). In this form, the Lemma shows 
that we can, in fact, teleport any 2x2 unitary matrix transformation U. We 
take M — U, and take the othogonal basis provided by the Lemma. Then a 
2-qubit classical transmission from Alice to Bob will enable Bob to identify 
the measured state and he can rotate it back to U\il)). 
Note that for H = V®^ we can consider the matrices 

where a = («(!), • • • , a{n)) and P = • • • , are bit-strings of length 

n and Tq^q = I, Tq i = X, Ti o = Y, Ti i = Z are the modified Pauli matri- 
ces discussed above. Then just as in the above Lemma, if C/ is a unitary 
matrix defined on H, then the set of measurement states {Tq.,i3^\ for the ma- 
trices Ta^pU are an orthogonal basis for H ® H. Hence we can teleport the 
action of the arbitrary unitary matrix U from Alice to Bob, at the expense 
of a transmission of 2" classical bits of information. This means that, we 
can construct an arbitrary unitary transformation (hence an idealized quan- 
tum computer) almost entirely by using quantum measurments. This result 
should be compared with the results of [20], and [42], which we shall do in 
a paper subsequent to the present work. If Alice and Bob conicide as ob- 
servers, then there is no need to transmit the classical bits. The result of a 
given measurement is an instruction to perform one of a preselected collection 
of unitary transformations on the resulting state. 

There are a number of lines that we will follow in susequent papers related 
to the points made in this section. In particular, it is certainly of interest 
that one can partially teleport transformations M that are not unitary, at the 
cost of having only partial information beforehand of the success of any given 
measurement. In particular, this means that we could consider computing 
results such as traces or generalized traces of matrices that are not unitary. In 



42 



this way we could examine computations of knot and link invariants that are 
based on non-unitary solutions to the Yang-Baxter equation. All of this will 
be the subject of another paper. In the next section we turn to the subject 
of quantum computation of link invariants based on unitary solutions to the 
Yang-Baxter equation. 



7 Unitary Representations of the Braid Group 
and the Corresponding Quantum Comput- 
ers 

Many questions are raised by the formulation of a quantum computer asso- 
ciated with a given link diagram configured as preparation, unitary transfor- 
mation and measurement. Unitary solutions to the Yang-Baxter equation (or 
unitary representations of the Artin braid group) that also give link invari- 
ants are not so easy to come by. Here we give a unitary representation that 
computes the Jones polynomial for closures of 3-braids. This representation 
provides a test case for the corresponding quantum computation. We now 
analyse this representation by making explicit how the bracket polynomial 
[22, 23, 28, 25] is computed from it. 

The idea behind the construction of this representation depends upon the 
algebra generated by two single qubit density matrices (ket-bras). Let \v) 

and \w) be two qubits in V, a complex vector space of dimension two over 
the complex numbers. Let P — \v){v\ and Q — \w){w\ be the corresponding 
ket-bras. Note that 

p2 = |f pP, 

PQP = \{v\w)\''P, 
QPQ^\{v\w)\^Q. 

P and Q generate a representation of the Temperley-Lieb algebra [25] . One 
can adjust parameters to make a representation of the three-strand braid 
group in the form 

si I — > rP + si, 
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S2 I — ^ tQ + ul, 

where / is the identity mapping on V and r, s, t, -u are suitably chosen scalars. 
In the following we use this method to adjust such a representation so that it 
is unitary. Note that it is possible for the representation to be unitary even 
though its mathematical "parts" P and Q are not unitary. Note also that the 
resulting representation is made entirely from local unitary transformations, 
so that while there is measurement of topological entanglement, there is no 
quantum entanglement of any corresponding quantum states. 

The representation depends on two symmetric but non-unitary matrices 
Ul and U2 with 



d 




and 



Uo 



d-^ Vl - d-^ 
VI - d-'^ d - d-^ 



Note that Ul = dUi and C/| = dU^. Moreover, U1U2U1 = Ui and U2U1U2 = 
Ul. This is an example of a specific representation of the Temperley-Lieb 
algebra [22, 25]. The desired representation of the Artin braid group is given 
on the two braid generators for the three strand braid group by the equations: 



^{si) ^ AI + A-^Ui, 

^{S2) = AI + A-^U2. 
Here / denotes the 2x2 identity matrix. 

For any A with d — —A"^ — A~^ these formulas define a representation of the 
braid group. With A — e*^, we have d — —2cos{29). We find a specific range 
of angles |^| < vr/G and |^ — 7r| < 71/6 that give unitary representations of the 
three-strand braid group. Thus a specialization of a more general represention 
of the braid group gives rise to a continuum family of unitary representations 
of the braid group. 

Note that tr{Ui) = triJJi) = d while tr{UiU2) = tr{U2Ui) = 1. If 6 is 
any braid, let I{h) denote the sum of the exponents in the braid word that 



44 



expresses b. For b a three-strand braid, it follows that 

$(6) = A^wj + n(6) 

where / is the 2x2 identity matrix and 11(6) is a sum of products in the 
Temperley-Lieb algebra involving Ui and U2- Since the Temperley-Lieb al- 
gebra in this dimension is generated by I,Ui, U2, U1U2 and U2U1, it follows 
that the value of the bracket polynomial of the closure of the braid 6, denoted 
<b>, can be calculated directly from the trace of this representation, except 
for the part involving the identity matrix. The result is the equation 

< & >= A^^''^(f + tr{U{b)) 

where b denotes the standard braid closure of b, and the sharp brackets denote 
the bracket polynomial. From this we see at once that 

<b>^trmb))+A^^''\d^ -2). 

It follows from this calculation that the question of computing the bracket 
polynomial for the closure of the three-strand braid b is mathematically 
equivalent to the problem of computing the trace of the matrix $(6). To 
what extent can our quantum computer determine the trace of this matrix? 
We have seen just before this subsection that a quantum computation can 
determine the absolute value of the trace by repeated trials. This shows that 
a major portion of the Jones polynomial for three strand braids can be done 
by quantum computation. 

7.1 The Invariant based on R 

A second example is given by the invariant discussed in the previous section. 
In that case, we have the formula 

r{b) = tr{repn{b)) 

taken up to multiples of the square root of 2, and the matrix repn{b) is 
unitary for any braid b in an n-strand braid group for arbitrary positive 
integer n. This invariant can be construed as the trace of unitary matrix for 
a quantum computation. Since, as we have seen, knowledge of the invariant 
often depends upon knowing the global sign of the trace of repn{b), it is not 
enuough to just compute the absolute value of this trace. Nevertheless some 
topological information is available just from the absolute value. 
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8 Quantum Entanglement and Topological 
Entanglement 

The second question about unitary solutions to the Yang-Baxter equation is 
the matter of understanding their capabihties in entanghng quantum states. 
We use the criterion that 

(f) = a|00) +6|01) + c|10) + 

is entangled if and only ii ad — be ^ 0. This criterion is generalized to higher 
dimensional pure states in our papers [26, 27, 37, 38, 40] . 

In [25, 26, 27], we discovered families of unitary solutions to the Yang-Baxter 
equation that detect topological linking if and only if the gates corresponding 
to these solutions can entangle quantum states. 

Is there a deeper connection between topological entanglement and quan- 
tum entanglement? We believe that more exploration is called for before a 
definitive answer to this question can be formulated. We need more bridges 
between quantum topology and quantum computation. 

The matrix 

/ 1/^2 1/V2 \ 

1/V2 -1/72 

1/V2 1/V2 

^ -1/V2 1/V2 , 

is a unitary solution of the Yang-Baxter equation; and it is highly entangling 
for quantum states. It takes the standard basis for the tensor product of two 
single qubit spaces to the Bell basis. On the topological side, R generates a 
new and non-trivial invariant of knots and links. On the quantum side, R is 
a universal gate at the same level as CNOT^ as we showed in Theorems 2 
and 3. Thus i? is a good example of a transformation that can be examined 
fruitfully in both the quantum and the topological contexts. 
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8.1 Linking Numbers and the Matrix R' . 



The unitary R' matrix that we have considered in this paper gives rise to a 
non-trivial invariant of hnks. The discussion in this section summarizes our 
treatemcnt of this invariant in [26]. Here we discuss the invariant associated 
with the speciahzation of E! with so that 

" a " 

fi"= 

c ■ 
.000a. 

The invariant is calculated from a state summation associated with the 
matrix R' and can be shown to have the form 

Zk = 2(1 + (cVa')'"^''^) 

for two-comonent links K, where lk{K) denotes the linking number of the 
two components of K. We show that for this specialization of the R' matrix 
the operator R! entangles quantum states exactly when it can detect linking 
numbers in the topological conteoct. 



Here is a description of the state sum: Label each component of the 
diagram with either or 1. Take vertex weights of a or c for each local la- 
belling of a positive crossing as shown in Figure 14. For a negative crossing 
(obtained by interchanging over-crossing and under- crossing segments at a 
positive crossing) the corresponding labels are 1/a and 1/c (which are the 
complex conjugates of a and c repsectively, when a and c are unit complex 
numbers). Let each state (labeling of the diagram by zeroes and ones) con- 
tribute the product of its vertex weights. Let S(i^) denote the sum over all 
the states of the products of the vertex weights. Then one can verify that 
Z{K) = a~'^^^^'E{K) where w{K) is the sum of the crossing signs of the 
diagram K. 
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Figure 14 - Positive Crossing Weights 



Figure 15 - Zero-One States for the Hopf Link 



For example, view Figure 15. Here we show the zero-one states for the 
Hopf hnk H. The 00 and 11 states each contribute a^, while the 01 and 10 
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states contribute Hence = 2(a^ + c^) and 

a-"'(^)S(i7) = 2(1 + {^lc?f) = 2(1 + (cVa')'^(^)), 

as expected. 

The calculation of the invariant in this form is an analysis of quantum 
networks with cycles in the underlying graph. In this form of calculation we 
are concerned with those states of the network that correspond to labelings by 
qubits that are compatible with the entire network structure. One considers 
only quantum states that are compatible with the interconnectedness of the 
network as a whole. 



8.2 The Question About Invariants and Entanglement 

We have seen that there are examples, such as the one given above, where 
topological entanglement measures, and measures of quantum entanglement 
are related to one another. In that example we found the the solution R' to 
the Yang-Baxter equation would, as an operator on states, entangle quantum 
states exactly when the invariant could measure linking numbers. We have 
also discussed the invariant assoicated with the universal gate it! and shown 
that it detects many topological situations that are quite subtle. For example, 
it can measure the linkedness of the Borromcan rings and the linkedness of 
the Whitehead link, both of which are situations where the linking numbers 
are zero. And yet, we have also given an example, in the previous section, of a 
representation of the braid group on three strands, B^, (not constructed from 
a solution to the Yang-Baxter equation) that produces the Jones polynomial 
for closures of three-stranded braids, but is defined on a single qubit. Since 
this last representation acts only on one qubit, there is no entanglement 
associated with it. Therefore it remains, at this writing, unclear just what is 
the relationship between the quantum entangling properties of braid group 
representations and their ability to measure topological entanglement. In 
a sequel to this paper we will concentrate this analysis just on invariants 
assoicated with solutions to the Yang-Baxter equation. 
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8.3 The Aravind Hypothesis 

Link diagrams can be used as graphical devices and holders of information. 
In this vein Aravind [1] proposed that the entanglement of a link should 
correspond to the entanglement of a state. Observation of a link would be 
modeled by deleting one component of the link. A key example is the Bor- 
romean rings. See Figures 7 and 16. 



Figure 16 - Borromean Rings 

Deleting any component of the Boromean rings yields a remaining pair of 
unlinked rings. The Borromean rings are entangled, but any two of them are 
unentangled. In this sense the Borromean rings are analogous to the GHZ 
state \GHZ) = (l/\/2)(|000) + |111)). Observation in any factor of the GHZ 
yields an unentangled state. Aravind points out that this property is basis 
dependent. We point out that there are states whose entanglement after an 
observation is a matter of probability (via quantum amplitudes). Consider 
for example the state 

IV') = (1/2)(|000) + |001) + |101) + |110)). 



Observation in any coordinate yields an entangled or an unentangled state 
with equal probability. For example 

|V^) = (1/2)(|0)(|00) + |01)) + |1)(|01) + |10))) 



so that projecting to |0) in the first coordinate yields an unentangled state, 
while projecting to |1) yields an entangled state, each with equal probability. 
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New ways to use link diagrams must be invented to map the properties 
of such states. We take seriously the problem of classifying the topological 
entanglement patterns of quantum states. We are convinced that such a clas- 
sification will be of practical importance to quantum computing, distributed 
quantum computing and relations with quantum information protocols. 



9 Braiding and Topological Quantum Field 
Theory 

The purpose of this section is to discuss in a very general way how braiding is 
related to topological quantum field theory and to the enterprise [17] of using 
this sort of theory as a model for anyonic quantum computation. The ideas in 
the subject of topological quantum field theory are well expressed in the book 
[2] by Michael Atiyah and the paper [50] by Edward Witten. The simplest 
case of this idea is C. N. Yang's original interpretation of the Yang-Baxter 
Equation [52] . Yang articulated a quantum field theory in one dimension of 
space and one dimension of time in which the i?-matrix (meaning here any 
matrix satisfying the Yang-Baxter equation) was regarded as giving the scat- 
tering ampitudes for an interaction of two particles whose (let us say) spins 
corresponded to the matrix indices so that is the amplitude for particles 
of spin a and spin b to interact and produce particles of spin c and d. Since 
these interactions are between particles in a line, one takes the convention 
that the particle with spin a is to the left of the particle with spin b, and the 
particle with spin c is to the left of the particle with spin d. If one follows the 
braiding diagram for a concatenation of such interactions, then there is an 
underlying permutation that is obtained by following the braid strands from 
the bottom to the top of the diagram (thinking of time as moving up the 
page). Yang designed the Yang-Baxter equation so that the amplitudes for 
a composite process depend only on the underlying permutation correspond- 
ing to the process and not on the individual sequences of interactions. The 
simplest example of this is the diagram for the Yang-Baxter equation itself 
as we have shown it in Figure 1. 

In taking over the Yang-Baxter equation for topological purposes, we 
can use the same intepretation, but think of the diagrams with their under- 
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and over-crossings as modeling events in a spacetime with two dimensions 
of space and one dimension of time. The extra spatial dimension is taken 
in displacing the woven strands perpendicular to the page, and allows us to 
use both braiding operators R and as scattering matrices. Taking this 
picture to heart, one can add other particle properties to the idealized the- 
ory. In particular one can add fusion and creation vertices where in fusion 
two particles interact to become a single particle and in creation one par- 
ticle changes (decays) into two particles. Matrix elements corresponding to 
trivalent vertices can represent these interactions. See Figure 17. 



Once one introduces trivalent vertices for fusion and creation, there is 
the question how these interactions will behave in respect to the braiding 
operators. There will be a matrix expression for the compositions of braid- 
ing and fusion or creation as indicated in Figure 18. Here we will restrict 
ourselves to showing the diagrammatics with the intent of giving the reader 
a flavor of these structures. It is natural to assume that braiding intertwines 
with creation as shown in Figure 19 (similarly with fusion). This intertwin- 
ing identity is clearly the sort of thing that a topologist will love, since it 
indicates that the diagrams can be interpreted as embeddings of graphs in 
three-dimensional space. Thus the intertwining identity is an assumption like 
the Yang-Baxter equation itself, that simplifies the mathematical structure 
of the model. 




Figure 17 -Creation and Fusion 
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Figure 18 - Braiding 




Figure 19 - Intertwining 

It is to be expected that there will be an operator that expresses the 
recoupling of vertex interactions as shown in Figure 20 and labeled by Q. 
The actual formalism of such an operator will parallel the mathematics of 
recoupling for angular momentum. See for example [24]. If one just considers 
the abstract structure of recoupling then one sees that for trees with four 
branches (each with a single root) there is a cycle of length five as shown in 
Figure 21. One can start with any pattern of three vertex interactions and 
go through a sequence of five recouplings that bring one back to the same 
tree from which one started. It is a natural simplifying axiom to assume that 
this composition is the identity mapping. This axiom is called the pentagon 
identity. 
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Figure 20 - Recoupling 




Figure 21 - Pentagon Identity 



Finally there is a hexagonal cycle of interactions between braiding, recou- 
pling and the intertwining identity as shown in Figure 22. One says that the 
interactions satisfy the hexagon identity if this composition is the identity 
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Figure 22 - Hexagon Identity 

A three-dimensional topological quantum field theory is an algebra of inter- 
actions that satisfies the Yang-Baxter equation, the intertwining identity, the 
pentagon identity and the hexagon identity. There is not room in this sum- 
mary to detail the remarkable way that these properties fit into the topology 
of knots and three-dimensional manifolds. As the reader can see, a three di- 
mensional TQFT is a highly simplified theory of point particle interactions 
in 2 -|- 1 dimensional spacetime. It can be used to articulate invariants of 
knots and links and invariants of three manifolds. The reader interested in 
the SU (2) case of this structure and its implications for invariants of knots 
and three manifolds can consult [24, 28, 32, 9, 43]. One expects that physical 
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situations involving 2 + 1 spacetime will be approximated by such an idealized 
theory. It is thought for example, that aspects of the quantum Hall effect 
will be related to topological quantum field theory [49]. One can imagine 
a physics where the space is two dimensional and the braiding of particles 
corresponds to their exchanges as though circulating around one another in 
the plane. Such particles that, unlike fcrmions, do not just change the am- 
plitude by a sign under interchange, but rather by a complex phase or even 
a linear combination of states are called anyons. It is hoped that TQFT 
models will describe applicable physics. One can think about the possible 
applications of anyons to quantum computing. The TQFT's then provide 
a class of anyonic models where the braiding is essential to the physics and 
to the quantum computation. We have given a sketch of this approach here 
to give the reader a picture of one of the possibilities of using braiding in 
quantum computing. 




Figure 23 - A More Complex Braiding Operator 
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The key point in the apphcation and relationship of TQFT and quantum 
information theory is, in our opinion, contained in the structure illustrated 
in Figure 23. There we show a more complex braiding operator, based on the 
composition of recoupling with the elementary braiding at a vertex. (This 
structure is implicit in the Hexagon identity of Figure 22.) The new braiding 
operator is a source of unitary representations of braid group in situations 
(which exist) where the recoupling transformations are themselves unitary. 
This kind of pattern is implicitly utilized in the work of Freedman and col- 
laborators [14, 15, 16, 17, 18] and in the case of classical angular momentum 
formalism has been dubbed a "spin-network quantum simlator" by Rasetti 
and collaborators [47]. 



10 Discussion 

It is natural to expect relationships between topology and quantum mechan- 
ics. For example, Dirac [11] described the relationship between an observer 
and a fermion by using the properties of twisted belts embedded in three 
dimensional space. These properties vividly portray the consequences of the 
fact that SU{2) double covers SO{S). The rotation group SO{S) and the 
unitary group SU (2) are involved since a rotation of the observer is mapped 
to a unitary transformation of the wave function. The topology of the belt 
gives a direct way to image the properties of this connection, with one full 
rotation changing the sign of the wave function, while two full rotations do 
not change that sign. In the topological picuture, that relationship between 
one object and another object rotated relative to the first object is depicted 
by a belt connecting them. Topological properties of the belt mimic the 
orientation - entanglement relation. 

How might such relationships between topology and quantum mechan- 
ics impinge upon quantum computing? The Dirac string trick suggests that 
topology may enter in the structure of non-locality and entanglement. On the 
quantum computing side, we know many uses for entangled states (e.g. tele- 
portation protocols); and one wants to understand the role of entanglement 
in the efficiency of computing procedures. Entanglement in quantum me- 
chanics and entanglement (linking and knotting) in topology can be related 
in a number of ways that give rise to a host of research questions. 
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We would like to state some general properties of this quest for relation- 
ship between topology and quantum mechanics: It is normally assumed that 
one is given the background space over which quantum mechanics appears. 
In fact, it is the already given nature of this space that can make non- locality 
appear mysterious. In writing |0) = (|01) + |10))/v^, we indicate the entan- 
gled nature of this quantum state without giving any hint about the spatial 
separation of the qubits that generate the first and second factors of the ten- 
sor product for the state. This split between the properties of the background 
space and properties of the quantum states is an artifact of the rarefied form 
given to the algebraic description of states, but it also points out that it is 
exactly the separation properties of the topology on the background space 
that arc implicated in a discussion of non-locality. 

Einstein, Podolsky and Rosen might have argued that if two points in 
space are separated by disjoint open sets containing them, then they should 
behave as though physically independent. Such a postulate of locality is 
really a postulate about the relationship of quantum mechanics to the topol- 
ogy of the background space. The Dirac string trick can be understood in 
a similar manner. In this way, we see that discussions of non-locality in 
quantum mechanics are in fact discussions of the relationship between prop- 
erties of the quantum states and properties of the topology of the background 
space. Subtle questions related to metric and change of metric give rise to 
the well-known problems of quantum gravity (since general relativity must 
take into account the subtleties of the spacetime metric and the topology of 
spacetime) . 

Approaches such as Roger Penrose's spin networks and the more recent 
work of John Baez, John Barrett, Louis Crane Lee Smolin, and others suggest 
that spacetime structure should emerge from networks of quantum interac- 
tions occurring in a pregeometric or process phase of physicality. In such a 
spin network model, there would be no separation between topological prop- 
erties and quantum properties. We intend to carry this discussion to the 
spin network or to the spin foam level. It is our aim to deepen the discussion 
of topology and quantum computing to a level where this can be done in a 
uniform manner. 

The spin network level is already active in topological models such as the 
Jones polynomial, the so called quantum invariants of knots, links and three- 
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manifolds, topological quantum field theories [2, 50], and related anyonic 
models for quantum computing [14, 15, 16]. For example, the bracket model 
[22, 23, 28, 25] for the Jones polynomial can be reahzed by generalization of 
the Penrose SU{2) spin nets to the quantum group SU{2)q. 

Since the advent of knot invariants such as the Jones polynomial, spin 
network studies have involved g-deformations of classical spin networks and 
the corresponding topological properties. These g-deformations are, in turn, 
directly related to properties of g-deformed Lie algebras (quantum groups, 
Hopf algebras) containing solutions to the Yang-Baxter equation. Solutions 
to the Yang-Baxter equation are maps R : V ®V — > V (^V on the tensor 
product of two vector spaces that represent topological braiding. 

A direct question important for us is the determination of unitary so- 
lutions to the Yang-Baxter equation, and the investigation of both their 
topological properties and their quantum information properties. For the 
latter we want to know what role such solutions (matrices) can play in quan- 
tum computing. Specific questions are how such a matrix can be used in a 
quantum computational model for a link invariant, and can the matrix in 
question map unentangled states to entangled states. Some of these specific 
phenomena have been discussed in this paper. 
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